REDHORSE-REcombination and Double crossover detection in Haploid Organisms using next-geneRation SEquencing data

BACKGROUND: Next-generation sequencing technology provides a means to study genetic exchange at a higher resolution than was possible using earlier technologies. However, this improvement presents challenges as the alignments of next generation sequence data to a reference genome cannot be directly used as input to existing detection algorithms, which instead typically use multiple sequence alignments as input. We therefore designed a software suite called REDHORSE that uses genomic alignments, extracts genetic markers, and generates multiple sequence alignments that can be used as input to existing recombination detection algorithms. In addition, REDHORSE implements a custom recombination detection algorithm that makes use of sequence information and genomic positions to accurately detect crossovers. REDHORSE is a portable and platform independent suite that provides efficient analysis of genetic crosses based on Next-generation sequencing data. RESULTS: We demonstrated the utility of REDHORSE using simulated data and real Next-generation sequencing data. The simulated dataset mimicked recombination between two known haploid parental strains and allowed comparison of detected break points against known true break points to assess performance of recombination detection algorithms. A newly generated NGS dataset from a genetic cross of Toxoplasma gondii allowed us to demonstrate our pipeline. REDHORSE successfully extracted the relevant genetic markers and was able to transform the read alignments from NGS to the genome to generate multiple sequence alignments. Recombination detection algorithm in REDHORSE was able to detect conventional crossovers and double crossovers typically associated with gene conversions whilst filtering out artifacts that might have been introduced during sequencing or alignment. REDHORSE outperformed other commonly used recombination detection algorithms in finding conventional crossovers. In addition, REDHORSE was the only algorithm that was able to detect double crossovers. CONCLUSION: REDHORSE is an efficient analytical pipeline that serves as a bridge between genomic alignments and existing recombination detection algorithms. Moreover, REDHORSE is equipped with a recombination detection algorithm specifically designed for Next-generation sequencing data. REDHORSE is portable, platform independent Java based utility that provides efficient analysis of genetic crosses based on Next-generation sequencing data. REDHORSE is available at http://redhorse.sourceforge.net/. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-015-1309-7) contains supplementary material, which is available to authorized users

Genetic mapping reveals that sinefungin resistance in Toxoplasma gondii is controlled by a putative amino acid transporter locus that can be used as a negative selectable marker

Quantitative trait locus (QTL) mapping studies have been integral in identifying and understanding virulence mechanisms in the parasite Toxoplasma gondii. In this study, we interrogated a different phenotype by mapping sinefungin (SNF) resistance in the genetic cross between type 2 ME49-FUDR(r) and type 10 VAND-SNF(r). The genetic map of this cross was generated by whole-genome sequencing of the progeny and subsequent identification of single nucleotide polymorphisms (SNPs) inherited from the parents. Based on this high-density genetic map, we were able to pinpoint the sinefungin resistance phenotype to one significant locus on chromosome IX. Within this locus, a single nonsynonymous SNP (nsSNP) resulting in an early stop codon in the TGVAND_290860 gene was identified, occurring only in the sinefungin-resistant progeny. Using CRISPR/CAS9, we were able to confirm that targeted disruption of TGVAND_290860 renders parasites sinefungin resistant. Because disruption of the SNR1 gene confers resistance, we also show that it can be used as a negative selectable marker to insert either a positive drug selection cassette or a heterologous reporter. These data demonstrate the power of combining classical genetic mapping, whole-genome sequencing, and CRISPR-mediated gene disruption for combined forward and reverse genetic strategies in T. gondii

Gene expression in Leishmania is regulated predominantly by gene dosage

ABSTRACT Leishmania tropica, a unicellular eukaryotic parasite present in North and East Africa, the Middle East, and the Indian subcontinent, has been linked to large outbreaks of cutaneous leishmaniasis in displaced populations in Iraq, Jordan, and Syria. Here, we report the genome sequence of this pathogen and 7,863 identified protein-coding genes, and we show that the majority of clinical isolates possess high levels of allelic diversity, genetic admixture, heterozygosity, and extensive aneuploidy. By utilizing paired genome-wide high-throughput DNA sequencing (DNA-seq) with RNA-seq, we found that gene dosage, at the level of individual genes or chromosomal “somy” (a general term covering disomy, trisomy, tetrasomy, etc.), accounted for greater than 85% of total gene expression variation in genes with a 2-fold or greater change in expression. High gene copy number variation (CNV) among membrane-bound transporters, a class of proteins previously implicated in drug resistance, was found for the most highly differentially expressed genes. Our results suggest that gene dosage is an adaptive trait that confers phenotypic plasticity among natural Leishmania populations by rapid down- or upregulation of transporter proteins to limit the effects of environmental stresses, such as drug selection. IMPORTANCE Leishmania is a genus of unicellular eukaryotic parasites that is responsible for a spectrum of human diseases that range from cutaneous leishmaniasis (CL) and mucocutaneous leishmaniasis (MCL) to life-threatening visceral leishmaniasis (VL). Developmental and strain-specific gene expression is largely thought to be due to mRNA message stability or posttranscriptional regulatory networks for this species, whose genome is organized into polycistronic gene clusters in the absence of promoter-mediated regulation of transcription initiation of nuclear genes. Genetic hybridization has been demonstrated to yield dramatic structural genomic variation, but whether such changes in gene dosage impact gene expression has not been formally investigated. Here we show that the predominant mechanism determining transcript abundance differences (>85%) in Leishmania tropica is that of gene dosage at the level of individual genes or chromosomal somy

Rhoptry proteins ROP5 and ROP18 are major murine virulence factors in genetically divergent South American strains of Toxoplasma gondii

Toxoplasma gondii has evolved a number of strategies to evade immune responses in its many hosts. Previous genetic mapping of crosses between clonal type 1, 2, and 3 strains of T. gondii, which are prevalent in Europe and North America, identified two rhoptry proteins, ROP5 and ROP18, that function together to block innate immune mechanisms activated by interferon gamma (IFNg) in murine hosts. However, the contribution of these and other virulence factors in more genetically divergent South American strains is unknown. Here we utilized a cross between the intermediately virulent North American type 2 ME49 strain and the highly virulent South American type 10 VAND strain to map the genetic basis for differences in virulence in the mouse. Quantitative trait locus (QTL) analysis of this new cross identified one peak that spanned the ROP5 locus on chromosome XII. CRISPR-Cas9 mediated deletion of all copies of ROP5 in the VAND strain rendered it avirulent and complementation confirmed that ROP5 is the major virulence factor accounting for differences between type 2 and type 10 strains. To extend these observations to other virulent South American strains representing distinct genetic populations, we knocked out ROP5 in type 8 TgCtBr5 and type 4 TgCtBr18 strains, resulting in complete loss of virulence in both backgrounds. Consistent with this, polymorphisms that show strong signatures of positive selection in ROP5 were shown to correspond to regions known to interface with host immunity factors. Because ROP5 and ROP18 function together to resist innate immune mechanisms, and a significant interaction between them was identified in a two-locus scan, we also assessed the role of ROP18 in the virulence of South American strains. Deletion of ROP18 in South American type 4, 8, and 10 strains resulted in complete attenuation in contrast to a partial loss of virulence seen for ROP18 knockouts in previously described type 1 parasites. These data show that ROP5 and ROP18 are conserved virulence factors in genetically diverse strains from North and South America, suggesting they evolved to resist innate immune defenses in ancestral T. gondii strains, and they have subsequently diversified under positive selection

Genetic Divergence of Toxoplasma gondii Strains Associated with Ocular Toxoplasmosis, Brazil

Brazilian strains of T. gondii differ from lineages in North America and Europe; these differences may underlie severe ocular disease

Selection at a single locus leads to widespread expansion of toxoplasma gondii lineages that are virulent in mice

The determinants of virulence are rarely defined for eukaryotic parasites such as T. gondii, a widespread parasite of mammals that also infects humans, sometimes with serious consequences. Recent laboratory studies have established that variation in a single secreted protein, a serine/threonine kinase known as ROPO18, controls whether or not mice survive infection. Here, we establish the extent and nature of variation in ROP18among a collection of parasite strains from geographically diverse regions. Compared to other genes, ROP18 showed extremely high levels of diversification and changes in expression level, which correlated with severity of infection in mice. Comparison with an out-group demonstrated that changes in the upstream region that regulates expression of ROP18 led to an historical increase in the expression and exposed the protein to diversifying selective pressure. Surprisingly, only three atypically distinct protein variants exist despite marked genetic divergence elsewhere in the genome. These three forms of ROP18 are likely adaptations for different niches in nature, and they confer markedly different virulence to mice. The widespread distribution of a single mouse-virulent allele among geographically and genetically disparate parasites may have consequences for transmission and disease in other hosts, including humans

Single Crystal Functional Oxides on Silicon

Single crystalline thin films of complex oxides show a rich variety of functional properties such as ferroelectricity, piezoelectricity, ferro and antiferromagnetism etc. that have the potential for completely new electronic applications (1-2). Direct synthesis of such oxides on Si remains challenging due to the fundamental crystal chemistry and mechanical incompatibility of dissimilar interfaces (3-16). Here we report integration of thin (down to 1 unit cell) single crystalline, complex oxide films onto Si substrates, by epitaxial transfer at room temperature. In a field effect transistor using a transferred Pb0.2Zr0.8TiO3 (PZT) layer as the gate insulator, we demonstrate direct reversible control of the semiconductor channel charge with polarization state. These results represent the realization of long pursued but yet to be demonstrated single crystal functional oxides on-demand on silicon